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DENSITY PDFS OF DIFFUSE GAS IN THE MILKY WAY 

Elly M. Berkhuijsen 1 and Andrew Fletcher 2 



Abstract. The probability distribution functions (PDFs) of the average 
densities of the diffuse ionized gas (DIG) and the diffuse atomic gas are 
close to lognormal, especially when lines of sight at \b\ < 5° and \b\ > 5° 
are considered separately. Our results provide strong support for the 
existence of a lognormal density PDF in the diffuse ISM, consistent 
with a turbulent origin of density structure in the diffuse gas. 



1 Introduction 

Simulations of the interstellar medium (ISM) have shown that, if isothermal tur- 
bulence is shaping the structure of the medium, the density distribution becomes 



lognormal ( Elmegreen & Scalo 20041 de Avillez & Breitschwerdt 2005 Wada & 



Norman 2007). The latter authors found that for a large enough volume, and for 



a long enough simulation run, the physical processes causing the density variations 
in the ISM in a galactic disc can be regarded as random and independent events. 
Therefore, the PDF of log(density) becomes Gaussian and the density PDF log- 
normal. The medium is inhomogeneous on a local scale, but in a quasi-steady 
state on a global scale. 

Little observational evidence exists to test the results of the simulations. To 
date lognormal PDFs have been derived for the HI column densities in the LMC 



(Wada et al. 2000), for emission measures perpendicular to the Galactic plane 
of the DIG in the Milky Way dHill et all 120071 120081), and for the local volume 



density of dust near stars within 400 pc of the Sun ( Gaustad & Van Buren 1993 ) 

X 



Here we present the PDFs of average volume densities of the DIG and the diffuse 
atomic gas in the Milky Way, and show that they are consistent with lognormal 
distributions as well. 

We investigated the density PDFs of the DIG and of diffuse atomic hydrogen 
gas (HI) in the solar neighbourhood, using the dispersion measures of 34 pulsars 



at known distances (Berkhuijsen & Miiller 20081 and HI absorption lines towards 
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Fig. 1. (a) Six panels showing PDFs of HI. Left panels: PDFs of all diffuse HI, right 
panels: PDFs of only warm HI. (b) Three panels showing PDFs at high latitude (|6| > 
5°) for two HI phases and the diffuse ionized gas. 



375 stars (Diplas & Savage, 1994a b). Further details and an extensive discussion 
of some of the results presented here are given in Berkhuijsen & Fletcher (2008). 



2 Results and discussion 

In Fig. [T]ja)(left panels) we present the PDFs of the average volume density of HI, 
( tihi ), for the full sample of 375 stars (top) and for the stars at low (\b\ < 5°) and 
high (\b\ > 5°) latitudes (middle and bottom panels) (Diplas & Savage 1994a). 
The PDF for all stars is approximately lognormal above log( nni) = —1, but there 
is a clear excess at lower densities reflected in the large reduced-^ 2 statistic (see 
Table[T]). The PDFs for low and high |6| have a lognormal shape but their maxima 
are shifted with respect to each other. The latter sample clearly causes the low- 
density excess in Fig.[lJa)(top left). This is also the case for the PDFs for the lines 
of sight (LOS) towards stars probing warm HI (Diplas & Savage 1994b) shown 



in Fig. 1(a) (right panels). The shifts between the maxima show the difference in 
mean density between low and high latitudes (see Table [l]). Note that all densities 
of the warm sample have (n) < 0.3 cm~ 3 . 

The dispersions of the PDFs of the high |6| samples are twice those of the low 
|6| samples. A possible explanation is that the low \b\ LOS cross more turbulent 
cells than at high |6| ( Vazquez-Semadeni fc Garcia| 2001). 



Table [T] shows that the dispersions of the PDFs of the warm HI are slightly 
smaller than for the full sample. The combination of warm and cool (denser) gas in 
the full sample increases the dispersion because the density range becomes larger. 
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Table 1. Lognormal fits to the PDFs of (tihi) and (n c ). The fitted function is 
Y — (\/27rcr) _1 exp[— (log 10 (n ) — fj,) 2 /2a 2 ], where Y is the normalised number of pul- 
sars or stars per bin, and the reduced chi-square value for each fit is given. 







Position of 


maximum 


Dispersion 




Area 


N 




(n) 


a 


x 2 


All stars 


375 


-0.57 ±0.03 


0.27 ±0.02 


0.34 ±0.03 


6.0 


\b\ < 5° 


213 


-0.52 ±0.02 


0.30 ±0.02 


0.24 ±0.02 


2.7 


\b\ > 5° 


162 


-0.80 ±0.02 


0.16 ±0.01 


0.50 ±0.02 


0.4 


Cool HI \b\ > 5° 


G4 


-0.40 ±0.05 


0.40 ±0.04 


0.27 ±0.03 


1.3 


Warm HI 


140 


-0.98 ±0.03 


0.10 ±0.01 


0.28 ±0.03 


2.5 


|6| < 5° 


42 


-0.89 ±0.01 


0.13 ±0.01 


0.18 ±0.01 


0.1 


\b\ > 5° 


98 


-1.07 ±0.03 


0.09 ±0.01 


0.33 ±0.03 


1.0 


DIG |b| > 5° 


34 


-1.83 ±0.02 


0.015 ±0.001 


0.22 ±0.02 


0.3 



This is clearly visible in Fig. [TJb): addition of the samples at high \b\ of cool HI 
(top) and warm HI (middle) gives the distribution for all high \b\ LOS shown in 
the bottom-left panel of Fig. [l](a) . However, the mixture of cool and warm HI only 
partly explains the difference in dispersion of the full samples at low and high |6|, 
because this difference also exists for the samples of warm HI. 

In the bottom panel of Fig. [ljb) we show the PDF of DIG, which is also 
lognormal, below that of the warm HI at the same The temperatures of the 
two components are similar and if the ionized and atomic gas are well mixed, one 
would expect their dispersions to be the same. However, the dispersion of the DIG 
sample is about 30 per cent smaller than that of the warm HI sample (see Table [T]). 
A plausible explanation for the difference, which is also consistent with the higher 
density of the maximum in the diffuse («hi) PDF, is that low density regions 
are more readily ionized than higher density gas and that the average degree of 
ionization of the diffuse gas is substantially lower than 50 per cent. We estimate 
the degree of ionization to be about 14 per cent, using the densities for (n e ) and 
( riH i ) in Table [l] Alternatively, the DIG could have a higher mean temperature 
than the warm HI but with a smaller temperature range; in the simulations of |de| 



Avillez & Breitschwerdt ( |2005[ ) high temperature gas indeed has a lower median 



density and smaller dispersion. 



3 Conclusions 



The results presented in Section [2] lead us to the following conclusions: 
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1. The density PDF of the diffuse ISM is approximately lognormal, but there 
is a significant low density excess in the best-fit distribution. 

2. The density PDFs of the diffuse ISM in the disk (|6| < 5°) and away from 
the disk > 5°) are lognormal, but the positions of their maxima and the 
dispersions differ. These differences produce the low density excess in the 
PDF for all lines of sight. 

3. Several effects seem to influence the shape of the PDF. An increase of the 
number of clouds/cells along the LOS causes a decrease in the dispersion, 
whereas a mixture of cool (dense) and warm (less dense) gas causes an in- 
crease in the dispersion. 

These results provide strong support for the existence of a lognormal density 
PDF in the diffuse (i.e. average densities of n < 1cm -3 ) ionized and neutral 
components of the ISM. In turn, the form of the PDFs is consistent with the 
small-scale structure of the diffuse ISM being controlled by turbulence. 
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